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Abstract 

We present measurements of branching fractions for four previously unobserved S-meson decays 
with an r/ or r]' meson in the final state. The data sample corresponds to 182 million BB pairs 
produced from e"'"e~ annihilation at the T(45) resonance. We measure the following branching 
fractions in units of 10"^: B{B^ ijuj) = 1.2 ± 0.6 ± 0.2 (< 2.3 at 90% C.L.), B{B^ r}K^) = 
2.5 ± 0.8 ± 0.1, B{B+ rjp+) = 8.6 ± 2.2 ± 1.1, and B{B+ r]'TT+) = 4.2 ± 1.0 ± 0.5, where the 
first error quoted is statistical and the second systematic. The charge asymmetries are Ach{B~^ 
'np+) = (7 ± 19 ± 2)% and Ach{B+ r/'7r+) = (24 ± 19 ± 1)%. AH resuhs are prehminary 
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1 INTRODUCTION 



We report results for measurements of the decay branching fractions of to the charmless final 
states^ rjio and r]K^, and of to rjp'^ and r/'yT"*". None of these decays have been observed 
definitively 0121 IB]- Measurements of the related decays — > i]K~^, B~^ — > r]TT~^ , and B r]'K 
were published recently [2111]. Charmless decays with kaons are usually expected to be dominated by 
b ^ s loop ("penguin") transitions, while h ^ u tree transitions are typically larger for the decays 
with pions and p mesons. However the B rjK decays are especially interesting since they are 
suppressed relative to the abundant B — > rj'K decays due to destructive interference between two 
penguin amplitudes [S]. The CKM-suppressed b ^ u amplitudes may interfere significantly with 
penguin amplitudes, possibly leading to large direct CP violation in B^ -qp^ and B^ — > r/'vr"'"|2]; 
numerical estimates are available in a few cases [71 |S1 . We search for such direct CP violation by 
measuring the charge asymmetry Ach = (r~ — r"'")/(r~ + F"*") in the rates = T{B^ — > /^), for 
each observed charged final state 

Charmless B decays are becoming useful to test the accuracy of theoretical predictions (HI HUl 
111! \V1\ \VA\ I14| I15| . Phenomenological fits to the branching fractions and charge asymmetries 
can be used to understand the importance of tree and penguin contributions and may provide 
sensitivity to the CKM angle 7 pT)] . 

2 THE BABAR DETECTOR AND DATASET 

The results presented here are based on data collected with the BABAR detector |17j at the PEP-II 
asymmetric e^e~ collider ^HI located at the Stanford Linear Accelerator Center. The samples 
come from an integrated luminosity of 166 fb~^ recorded at the T(45) resonance (center-of-mass 
energy ^/s = 10.58 GeV). This corresponds to 182 it 2 million BB pairs. 

Charged particles from the e^e~ interactions are detected and their momenta measured by a 
combination of a vertex tracker (SVT) consisting of five layers of double-sided silicon microstrip 
detectors and a 40-layer central drift chamber, both operating in the 1.5-T magnetic field of a super- 
conducting solenoid. We identify photons and electrons using a CsI(Tl) electromagnetic calorime- 
ter (EMC). Further charged particle identification (FID) is provided by the average energy loss 
(dE/dx) in the tracking devices and by an internally reflecting ring imaging Cherenkov detector 
(DIRC) covering the central region. 

3 SAMPLE SELECTION 

We reconstruct t], t]', uj, p'^ , p^, tt^, and K'^ candidates through their decays 7/ — > 77 (7777), 
r] 7r+7r"7r° (773^), rj' r/Tr+vr- (r/^^^), r/' p^j (r/^^), u 7r+7r"7r°, p+ 7r+7r°, p° tt+vt-, 
— > 77, and Kg vr'^vr". We make the requirements given in Tabled on the invariant mass of 
these particles' final states. For the r], oj, and rj' invariant masses these requirements are set loose 
enough to include sidebands, as these mass values are treated as observables in the maximum- 
likelihood (ML) fit described below. For K'^ candidates we further require the three-dimensional 
flight distance from the beam spot to be greater than three times its uncertainty in a fit that 
requires consistency between the flight and momentum directions. For modes with i? ^ — > 77 

^Except as noted explicitly, we use a particle name to denote either member of a charge conjugate pair. 
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Table 1: Selection requirements on the invariant mass of B daughter resonances (in MeV). 



State 


Requirement 




4yu < m(77j < buu 


Van 


OZU < Tn^TTTTTT ) < 0/U 


VrjiriT 


910 < mirjTTir) < 1000 


Vfry 


910 < m{pj) < 1000 


U 


735 < m(7r7r7r) < 825 


P+ 


470 < mivrvr) < 1070 


/ 


510 < m{TTn) < 1060 




120 < m(77) < 150 




486 < m(7r7r) < 510 



we impose a mode-dependent requirement on the decay angle to reject backgrounds reconstructed 
as very asymmetric decays. 

We make several PID requirements to ensure the identity of the charged pions and kaons. 
Secondary pions in ryajr, r/', and u candidates are rejected if their DIRC, dE/dx, and EMC outputs 
satisfy tight consistency with kaons, protons, or electrons. For the decays to an rj' meson and 
a charged pion or kaon, the latter (primary) track must have an associated DIRC signal with a 
Cherenkov angle within 3.5 standard deviations (cr) of the expected value for either the tt or K 
hypothesis. The discrimination between pion and kaon primary tracks is treated in the ML fit. 

The number of candidates found per event is at or below about 1.10 for all modes except for 
those with the final states ??37r/o and ??37rW where it is about 1.3. We choose the candidate whose 
daughter resonance mass(es) lie nearest the expected mean value. 

A i?-meson candidate is characterized kinematically by the energy-substituted mass ttt-es = 
[{^s + Po • pb)'^/Eq — p^] 2 and energy difference AE = E*g — \^/s, where the subscripts and B 
refer to the initial T(45) and to the B candidate, respectively, and the asterisk denotes the T{AS) 
frame. The resolution on IS.E (ttt-es) is about 30 MeV (3.0 MeV). We require |A£'| < 0.2 GeV and 
5.25 < rn-Es ^ 5.29 GeV, and include both of these observables in the ML fit. 

4 BACKGROUNDS 

Backgrounds arise primarily from random combinations in continuum e'^e~ qq events {q = 
u,d,s,c). We reject these by using the angle 9t in the T(45) frame between the thrust axis of 
the B candidate and that of the rest of the charged tracks and neutral calorimeter clusters in the 
event. The distribution of | cos^tI is sharply peaked near 1.0 for combinations drawn from jet-like 
qq pairs, and nearly uniform for i?-meson pairs. We require | cos 0t| < 0.9 for all modes except the 
high-background B^ v'pj'^^ decay. For this mode we determined that the sensitivity is maximal 
with I cos^tI < 0.65, based on the expected signal yield and its background-dominated statistical 
error. In the ML fit we also use a Fisher discriminant J-" |19j that combines four variables defined 
in the T{4S) frame: the angles with respect to the beam axis of the B momentum and B thrust 
axis, and the zeroth and second angular moments Lo,2 of the energy flow about the B thrust axis. 
The moments are defined by Lj = J2iPi ^ |cos0i|-' , where 9i is the angle with respect to the B 
thrust axis of track or neutral cluster i, pi is its momentum, and the sum excludes the B candidate 
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daughters. 

For the r/ ^ 77 modes we use additional event-selection criteria to reduce BB backgrounds 
from several charmless final states. We reduce background from B vr+vr'', K'^n^, and K^tt^ by 
rejecting rj^^ candidates that share a photon with any vr'^ candidate having momentum between 
1.9 and 3.1 GeV/c in the T(45) frame. Additionally, for B° rjK^ we require E* < 2.4 GeV to 
suppress background from B K*^ and related radiative-penguin decays. 

We use Monte Carlo (MC) simulation j2Uj for an initial estimate of the residual charmless 
BB background. Most of the contribution from b ^ c decays has a dependence on the ML fit 
observables that is similar to that for continuum events, and thus can be modeled as part of the 
continuum component. With a survey from MC we identify the few (mostly charmless) decays that 
may survive the candidate selection. We find these contributions to be negligible for several of our 
modes. Where they are not we include a component in the ML fit to account for them. 



5 MAXIMUM LIKELIHOOD FIT 



We obtain yields and Ach from extended unbinned maximum-likelihood fits with input observables 
AE, Jn-Es, and mres (the mass of the r], r]' , p'^, or to candidate). For the lo decays we also use 
TC = I cos 6h\, and for charged modes the PID variable St^^k- The helicity angle 6h is defined as 
the angle, measured in the uj rest frame, between the normal to the to decay plane and the flight 
direction of the oj with respect to its parent B. We incorporate PID information by using S-j^ (Sk), 
the number of standard deviations between the measured Cherenkov angle and that expected for 
pions (kaons). 

For each event i, hypothesis j (signal, continuum background, BB background), and flavor 
(primary Tr~^ or K^) k, we define the probability density function (PDF) 



Pj, = Vj{mEs')r,iAEl)V,ir)V,\ 



(1) 



The terms in brackets for S and Ti. pertain to modes with a primary charged track or lo daughters, 
respectively. The absence of correlations among observables in the background Pj^ is confirmed in 
the (background-dominated) data samples entering the fit. For the signal component, we correct 
for the effect of the neglect of small correlations (see below). The likelihood function is 



N 



jk 



(2) 



where Yjk is the yield of events of hypothesis j and flavor k found by maximizing £, and is the 
number of events in the sample. 

For the signal and BB background components we determine the PDF parameters from sim- 
ulation. For the continuum background we use {m-^,^^ AE) sideband data to obtain initial values, 
before applying the fit to data in the signal region, and ultimately by leaving them free in the final 
fit. We parameterize each of the functions 'Psig("iEs)) 'Psig(A£'fc), T'j{J^), Vj{Sk) and the peaking 
components of Vj{mj-cs) with either a Gaussian, the sum of two Gaussians or an asymmetric Gaus- 
sian function as required to describe the distribution. Slowly varying distributions (mass, energy or 
helicity-angle for combinatorial background) are represented by linear or quadratic dependencies. 
The peaking and combinatorial components of the u mass spectrum each have their own TC shapes. 
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The combinatorial background in /tt-es is described by the function x\/l — exp [—■^(1 — x^)] , with 
X = 2mEs/\/s and parameter Large control samples of B decays to charmed final states of sim- 
ilar topology are used to verify the simulated resolutions in AE and m-Es- Where the control data 
samples reveal differences from MC in mass or energy offset or resolution, we shift or scale the 
resolution function used in the likelihood fits. 

Before applying the fitting procedure to the data to extract the signal yields we subject it to 
several tests. Internal consistency is checked with fits to ensembles of "experiments" generated 
by Monte Carlo from the PDFs. Prom these we establish the number of parameters associated 
with the PDF shapes that can be left free in addition to the yields. Ensemble distributions of the 
fitted parameters verify that the generated values are reproduced with the expected resolution. 
The ensemble distribution of In C itself provides a reference to check the goodness of fit of the final 
measurement once it has been performed. 

We evaluate possible biases from our neglect of correlations among discriminating variables 
in the PDFs by fitting ensembles of simulated experiments into which we have embedded the 
expected number of signal events randomly extracted from the fully simulated MC samples. We 
find a positive bias of < 1 event for the tjuj and tjK^ modes. For rjp^ and 7]'tt~^ it is 7 to 17 events 
(~ 15% of the yield). Events from a weighted mixture of simulated BB background decays are 
included where significant, and so the bias we measure includes the effect of crossfeed from these 
modes. 



6 FIT RESULTS 

Free parameters of the fit include signal and background yields, background PDF parameters, and 
for charged modes the signal and background Ach- The free background PDF parameters are mean, 
width, and skewness for J^, slope for AE, slope of the combinatorial component and peak fraction 
for resonance mass, and ^ for ttt-es- 

The branching fraction for each decay chain is obtained from 

where Y is the yield of signal events from the fit, Yb is the fit bias discussed in the previous section, 
e is the efficiency, Yl is the product of daughter branching fractions that were forced to unity in 
the determination of e, and Nb is the number of produced B^ or B^ mesons. In Table HI we show 
for each decay mode the measured branching fraction together with the event yields and efficiencies. 
We assume that the decay rates of the T{AS) to B'^B~ and B^B^ are equal. The estimated purity 
is the ratio of the signal yield to the effective background plus signal; the sum of effective bkg plus 
signal is represented by the square of the uncertainty of the signal yield. 

In Figs. we show projections onto rn-ES and AE of subsamples enriched with a mode- 
dependent threshold requirement on the signal likelihood (computed without the PDF associated 
with the variable plotted). 



7 STATISTICAL AND SYSTEMATIC UNCERTAINTIES 

The statistical error on the signal yield and Ach is taken as the change in value that corresponds 
to an increase of — 21n£ by one unit from its minimum. The significance is taken as the square 
root of the difference between the value of — 21n£ (with systematic uncertainties included) for zero 
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Table 2: Signal yield Y, estimated purity P, detection efficiency e, daughter brancliing frac- 
tion product, significance S (with systematic uncertainties included), measured branching fraction, 
background {A^D and signal (Ach) charge asymmetries for each mode. 
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Figure 1: Projections of the B candidate rriEs (left) and Ai? (right) for — > t/u. Points with 
errors represent data, solid curves the full fit functions, and dashed curves the background functions. 
These plots are made with a requirement on the likelihood and thus do not show all events in the 
data sample. 

signal and the value at its minimum. For r/a; we quote a 90% confidence level (C.L.) upper limit, 
taken to be the branching fraction below which lies 90% of the total of the likelihood integral in 
the positive branching fraction region. For the charged modes we also give the charge asymmetry 

Ach- 

For the B^ — > r/'/i"*" fits we obtain yields also for the B^ — > r]'K^ decays. For both submodes 
these yields are consistent with the expectation from our previous measurements jl]. 
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25 5.26 5.27 5.28 5.29 

m^s (GeV) 




0.1 0.2 

AE (GeV) 



Figure 2: Projections of the B candidate m-Es (left) and /S.E (right) for rjK^. Points 

with errors represent data, sohd curves the full fit functions, and dashed curves the background 
functions. These plots are made with a requirement on the likelihood and thus do not show all 
events in the data sample. 




5.25 5.26 5.27 5.28 5.29 

(GeV) 




0.1 0.2 
AE (GeV) 



Figure 3: Projections of the B candidate mss (a) and AE (b) for B^ VP'^ ■ Points with errors 
represent data, solid curves the full fit functions, and dashed curves the background functions. 
These plots are made with a requirement on the likelihood and thus do not show all events in the 
data sample. 
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5.25 5.26 5.27 5.28 5.29 -0.2 -0.1 o!l 0.2 

nigs (GeV) AE (GeV) 

Figure 4: Projections of the B candidate mss (a) and AE (b) for y/vr"*". Points with 

errors represent data, sohd curves the full fit functions, and dashed curves the continuum plus 
— > rj'K~^ background functions. These plots are made with a requirement on the likelihood and 
thus do not show all events in the data sample. 



Most of the systematic uncertainties arising from lack of knowledge of the PDFs have been 
included in the statistical error since most background parameters are free in the fit. For the 
signal the uncertainties in PDF parameters are estimated from the consistency of fits to MC and 
data in control modes. Varying the signal PDF parameters within these errors, we estimate the 
uncertainties in the signal PDFs to be 1-8 events, depending on the mode. The uncertainty in the 
fit bias correction is taken to be half of the correction itself. Similarly we estimate the uncertainty 
from modeling the BB backgrounds by taking half of the contribution of that component to the 
fitted signal yield. These additive systematic errors are small for the rju and rjK^ modes, but 
dominant for ijp^ and r]TT^ . 

Uncertainties in our knowledge of the efficiency, found from auxiliary studies, include 0.8Nt%, 
l.bN^Vo, and 3.4% for a Kg decay, where Nt and are the number of signal tracks and photons, 
respectively. Our estimate of the B production systematic error is 1.1%. Published data pij 
provide the uncertainties in the i3-daughter product branching fractions (1%). The uncertainties 
in the efficiency from the event selection are 1% (3% in B'^ r]'p^Ti^) for the requirement on cos 0t 
and ~1% for PID. Using several large inclusive kaon and i?-decay samples, we find a systematic 
uncertainty for Ach of 1.1%, due mainly to the dependence of the reconstruction efficiency on the 
charge, for the high momentum pion from B^ rjir^ . The corresponding number for the softer 
charged pion from the p in B'^ r]p^ is 2%. The values of (see Table provide conffimation 
of this estimate. 

The pairs of separate daughter-decay measurements for each mode are combined by adding the 
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values of — 21n£ as functions of branching fraction, taking proper account of the correlated and 
uncorrelated systematic errors. We show these curves in Fig. [HI 

8 SUMMARY OF RESULTS 

In summary, we report preliminary results of searches for four charmless i?-meson decays. We find 
significant signals for the previously- undetected — > rjK^, — > r/p"*", and — > rl'K^ . The 
measured branching fractions are 

(1.2 ± 0.6 ± 0.2) X 10"^ (< 2.3x10"^), 
(2.5 ±0.8 ±0.1) X 10~^ 
(8.6 ±2.2 ± 1.1) X 10~^, 
(4.2 ±1.0 ±0.5) X 10"^, 

where the first error quoted is statistical, the second systematic; the upper limit is taken at 90% 
CL. For the B^ modes the charge asymmetries, are 

Aft(w+) = (7 ±19 ±2)%, 
-4e/.(V^+) = (24±19±1)%. 

Theoretical approaches to the study of these decays include those based on flavor SU(3) relations 
among many modes |9lllO | [T5 | ll6j. effective Hamiltonians with factorization and specific S-to- light- 
meson form factors |11| . perturbative QCD |12j . and QCD factorization [51 1131 [Hj . Our branching 
fraction measurements are generally in agreement with the ranges of these theoretical estimates. 
From global fits to the growing body of data on charmless B decays the component amplitudes and 
theoretically uncertain parameters of these models are coming to be significantly over-constrained 
[1() |I14| [TH]. Our measurement oi Ach in B^ — > r/'vr^ excludes the larger-magnitude negative values 
among the theoretical estimates |7| [H] . 
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B{R^ 7]io) = 

B{B^^ilK^) = 

B{B+ rip+) = 

B(B+ ri'Tr+) = 
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BF(B^^ Tj p^) X 10^ BF(B^^ r|'7i^) x 10^ 



Figure 5: Plots of individual and combined — 21n£ for branching fraction fits are shown for 
the decay rjuj (upper left), — > rjK^ (upper right), B^ — > 77/9+ (lower left), and B^ 

r/'vr"^ (lower right). Each plot shows the daughter modes as curves that are dashed {B^ — > rjy^Lo, 
B^ ^ r/3^K0, B+ ^ r/^^p+, B+ ^ ??;,,vr+) or dotted (5° ^ rys.w, 5° ^ r?-,^i^°, B+ ^ 7737^/?+, 
B^ rj' TT^), and the result of combining these as a solid curve. 
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